FORCES   ON  A  HORIZONTAL  CYLINDER 
DUE   TO   NON-LINEAR  WAVES 


Fred   Herman   Gehrman 


Library 

ite  School  Naval  Postgraduate  Sv. 

California  93940  Monterey,  California  93< 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 


FORCES  ON  A 

HORIZONTAL  CYLINDER 

DUE  TO 

NON-LINEAR  WAVES 

by 

Fred  He 

rman  Gehrman  Jr. 

Thes 

is  Advisor: 

C.J.  Garrison 

December  1972 


" 


kppfiovnd  {^oh.  pubtic.  n.<itQjOLi> e;   du>tAAbuutioyi  antiirUXzd. 


Forces  on  a  Horizontal  Cylinder  Due 
to  Non-Linear  Waves 


by 


Fred  Herman  ,Gehrman,  Jr. 
Lieutenant,  United  States  Navy 
B.  S.,  The  Creighton  University,  1966 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 


from  the 

NAVAL  POSTGRADUATE  SCHOOL 
December  1972 


ft* 


G  * 


Library 

Naval  Postgraduate  School 

Monterey,  California  93940 


ABSTRACT 

A  horizontal  cylinder,  located  near  the  floor  of  a  two 
dimensional  wave  channel  was  subjected  to  a  train  of  non- 
linear gravity  waves.   The  horizontal  and  vertical  forces 
were  measured  and  presented  in  dimensionless  form.   Experi- 
mental values  of  horizontal  and  vertical  force  coefficients 
are  presented  as  functions  of  dimensionless  wave  height  and 
dimensionless  wave  period.   The  dimensionless  force  coeffi- 
cients predicted  by  a  modified  Morrison  equation  and  a 
Froude-Krilov  force  are  compared  to  experimental  data. 
Fluid  particle  velocity  and  acceleration  values  were  calcu- 
lated from  Stokes  fifth-order  gravity  wave  theory.   Experi- 
mental dimensionless  wave  periods  from  10  to  200  were 
investigated. 

The  horizontal  force  coefficients  were  found  to  vary 
linearly  with  wave  height  for  dimensionless  period  values 
from  60  to  120.   The  vertical  force  coefficients  were 
found  to  be  inertia  dominated  at  low  dimensionless  wave 
heights  and  dominated  by  a  lift  force  at  higher  wave  heights. 
The  theory  predicted  experimental  values  of  force  coefficients 
with  good  accuracy,  especially  at  greater  water  depths. 
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I.   INTRODUCTION 

Exploration  for  and  use  of  ocean  resources  has  received 
increasing  attention  in  the  scientific  community  in  recent 
years.   A  considerable  effort  has  resulted  in  the  development 
of  technology  in  areas  such  as  off-shore  oil  exploration, 
recovery  of  petroleum  and  deployment  of  underwater  habitats. 
Oil  production  often  involved  the  use  of  very  large  submerged 
oil  storage  tanks  and  associated  piping  systems.   In  some 
locations,  oil  is  transferred  long  distances  in  submerged 
piping  systems,  often  leading  through  the  surf  zone  to 
facilities  ashore.   Also,  waste  disposal  systems  usually 
involve  deployment  of  large  diameter  outfall  pipes  which  are 
laid  through  the  surf  zone.   Proposals  for  large  scale  mining 
and  food  gathering  activities  indicate  the  trend  toward 
design  and  construction  of  large  submerged  structures. 

This  activity  has  generated  considerable  interest  in  the 
interaction  of  gravity  waves  with  submerged  objects  and 
particularly  in  the  forces  induced.   Most  previous  investi- 
gations [1,2,3]  in  this  area  center  on  the  use  of  a  so-called 
Morrison  Equation  [3],  which  involves  both  a  drag  and  iner- 
tial  component  of  force.   Unfortunately,  the  drag  and  added 
mass  coefficients  are  found  to  be  neither  constant  nor 
simply  shape  -  or  Reynolds  number  -  dependent  in  the  unsteady 
flow  that  occurs  in  the  interaction  of  waves  with  fixed 
bodies.   In  addition,  where  the  amplitude  of  wave  motion  is 
large  compared  to  the  body  dimensions,  flow  separation  occurs, 


and  the  coefficients  vary  as  the  flow  in  a  given  direction 
develops.   Finally,  for  small  depth  to  characteristic  dimen- 
sion ratios,  the  proximity  of  the  free  surface  may  influence 
the  value  of  the  coefficients,  even  if  no  flow  separation 
occurs. 

In  spite  of  these  difficulties,  the  Morrison  Equation 
will,  with  certain  restrictions,  provide  a  complete  descrip- 
tion of  the  horizontal  forces  acting  on  a  large  horizontal 
cylinder  in  contact  with  a  plane  boundary.   If  the  cylinder 
is  large  in  comparison  to  the  amplitude  of  the  fluid  motion, 
separation  does  not  occur,  and  the  resulting  wave-induced 
oscillatory  flow  may  be  considered  an  unseparated  potential 
flow  about  the  cylinder.   This  has  been  verified  by  Kuelegan 
and  Carpenter  [2]  and  by  Sarpkaya  and  Garrison  [4].   For 
such  cases,  the  drag  component  of  force  is  negligible,  while 
the  added  mass  coefficient  is  a  constant. 

Application  of  the  Morrison  Equation  to  determine  the 
vertical  components  of  force  for  a  bottom-mounted  cylinder 
yields  a  vertical  force  component  of  zero,  an  obviously 
invalid  result.   It  is  recognized  that  the  vertical  force 
should  be  composed  of  two  components;  one  associated  with 
the  inertia  effects  and  a  second  which  accounts  for  the  lift 
force  caused  by  the  increased  velocity  and  hence  decreased 
the  pressure  over  the  top  of  the  cylinder. 

The  inertial  component  of  vertical  force  can  be  approxi- 
mated by  what  is  commonly  known  as  a  Froude-Krilov  force. 
This  is  the  force  caused  by  the  pressure  distribution  around 
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the  surface  of  the  cylinder  existing  in  the  wave  if  the 
cylinder  was  not  present.   This  is  an  approximation  which 
does  not  account  for  the  presence  of  the  body  and  is 
therefore  considered  to  be  an  underestimate  of  the  force. 

The  lift  component  of  vertical  force  can  be  approximated 
by  use  of  the  unseparated,  potential  flow  model  used  in  the 
analysis  of  the  horizontal  forces,  i.e.,  that  of  a  uniform 
flow  past  a  cylinder  in  contact  with  a  rigid  boundary.   The 
lift  coefficient  for  this  case  has  been  determined  by  Dalton 
and  Helfenstein  [5]. 

A  study  carried  out  by  Johnson  [6]  for  horizontal  forces 
acting  on  a  bottom-mounted,  horizontal  cylinder  in  long  waves, 
assumed  that  at  long  wave  lengths,  the  horizontal  forces 
depended  on  wave  height  and  water  depth  only.   This  assumption 
is,  however,  of  doubtful  validity  in  the  range  of  Johnson's 
test. 

Shiller  [7]  measured  wave  forces  on  a  submerged  horizontal 
cylinder  due  to  small  amplitude  waves.   Over  the  range  of 
wave  heights  considered,  the  magnitude  of  the  horizontal 
force  was  found  to  vary  linearly  with  wave  height.   The 
magnitude  of  the  vertical  force  was  found  to  increase  in 
proportion  to  the  wave  height  squared. 

Perkinson  [8]  extended  Shiller' s  investigation  to  include 
larger  wave  heights  and  periods.   His  study  indicates  that 
horizontal  forces  increased  linearly  with  wave  heights  and 
becomes  independent  of  wave  period  for  large  wave  periods. 
The  vertical  forces  were  found  to  contain  two  regimes:  the 
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lower  wave  lengths  where  lift  force  predominates;  and  the 
higher  wave  lengths  where  inertial  force  is  dominant. 

In  order  to  evaluate  any  forces  using  the  above  analyses, 
it  is  necessary  to  first  evaluate  fluid  velocities,  accelera- 
tions and  pressure  by  use  of  some  wave  theory.  Fairly  simple 
expressions  may  be  derived  for  the  forces  by  use  of  Airy  wave 
theory.  However,  for  finite  amplitude  waves  in  finite  depth 
water,  a  non-linear  wave  theory  is  considered  to  be  necessary 
Stokes  fifth-order  wave  theory  has  been  chosen  for  this 
study. 

Stokes  originally  developed  a  second-order  theory  [9]  for 
the  case  of  a  non-linear  wave  in  water  of  finite  depth.   This 
method  has  been  extended  by  Borgman  and  Chappelear  [10]  to 
third  order.   Skjelbreia  and  Hendrickson  [11]  have  extended 
the  solution  to  third  and  fifth  order.   Bretschneider  [12] 
has  presented  a  method  for  extension  to  any  order. 

The  purpose  of  this  study  is  to  compare  the  analytical 
results  of  Stoke' s  fifth-order  wave  theory  with  experimental 
values  at  longer  wave  periods.  The  transition  section  of  the 
wave  channel  used  by  Perkinson  was  extended  approximately 
forty  feet  to  provide  a  slope  of  1:20.   This  provided  the 
possibility  of  extending  Perkinson' s  work  to  include  longer 
wave  periods.   This  study,  then,  is  an  extension  of  the 
work  of  Shiller  [7]  and  Perkinson  [8]. 
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II.   THEORETICAL  ANALYSIS 

A.  PROBLEM  DEFINITION 

The  problem  under  consideration  is  illustrated  in  Figure 
1.   A  train  of  gravity  waves  is  considered  to  progress  in 
the  positive  x-direction  in  water  of  depth  h,  with  rate  of 
propagation  C.   The  fluid  particle  velocity  is  expressed  in 
terms  of  the  horizontal  and  vertical  components,  u  and  v.   It 
is  of  primary  interest  to  determine  the  horizontal  and 
vertical  components  of  wave  force  acting  on  the  horizontal 
cylinder  which  is  in  contact  with  the  rigid  bottom. 

B.  DIMENSIONAL  ANALYSIS 

An  exact  analytical  solution  to  this  problem  is  quite 
formidable  and,  therefore,  an  approximate  solution  is  con- 
sidered in  this  work.   However,  first  it  is  instructive  to 
carry  out  a  dimensional  analysis  of  all  pertinent  parameters 
involved. 

In  general,  the  maximum  of  the  wave  force  per  unit 
length  acting  on  a  submerged  cylinder  due  to  wave  motion  is 
known  a  priori  to  depend  on  the  following  parameters: 

FmaxA  =  Vh'H'T'L'a>p'g>y)  (1) 


where 


F  =  wave  force 
h  =  water  depth 
H  =  wave  height 
T  =  wave  period 
L  =  wave  length 
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T=  Wave  Period 


//////////  /    /     /    /     /  7 — 7-7 — 7-7 


Fig.  I.    DEFINITION  OF  GEOMETRY 
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a  =  cylinder  radius 

p  =  fluid  density 

g  =  gravitational  constant 

u  =  fluid  viscosity 

%  =   cylinder  length 

However,  a  relationship  exists  between  the  parameters 

associated  with  the  incident  wave  (i.e.  h,H,T  and  L). 

i 

Consequently,  only  three  of  these  parameters  (e.g.  h,H  and 
T)  are  needed  to  completely  specify  the  incident  wave.   A 
dimensional  analysis  of  the  parameters  in  equation  (1)  yields 
the  following  groups: 

Pmax/pga2£=  f2(gT2/h,  h/a,  H/2a,  y/S,"Vgha2)   (2) 

The  last  term  on  the  right  side  of  equation  (2)  repre- 
sents the  ratio  of  Froude  number  to  Reynolds  number,  indi- 
cating the  ratio  of  viscous  to  inertial  forces.   It  is 
believed  that,  if  this  number  is  small,  it  may  be  neglected 
from  further  consideration,  and  equation  (2)  may  be  rewritten 

Vx  ■  VgT2/h>  5>  S)  (3) 

a.  p 

where  F  =  F        /pga   I 

max         max  Mfe 

H  =  H/2a 

h  =  h/a 

The  approximate  analytical  approach  of  this  study  seeks: 
1.   To  develop  expressions  for  the  fluid  particle  motion 
in  the  gravity  wave  as  functions  of  the  incident  wave 
parameters,  and 
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2).   To  develop  expressions  for  the  forces  on  the  sub- 
merged cylinder  as  functions  of  the  fluid  particle 
motion. 
The  assumptions  upon  which  this  analysis  is  based  are: 
1).   Water  will  be  considered  inviscid  and  incompressible, 

2).   The  cylinder  radius  will  be  considered  small  in 
comparison  to  the  dimensions  (h  and  L)  of  the 
incident  wave. 

C.   WAVE  THEORY 

The  motion  of  fluid  particles  in  a  gravity  wave  based  on 
assumptions  1)  and  2)  above  is  specified  by  the  following 
boundary  value  problem. 

1).   Governing  Equation 

Assumption  1)  implies  that  Laplace's  Equation 
governs  the  fluid  motion.   That  is: 

V2*  =  0  (4a) 

where  $  is  the  velocity  potential  and  the  two  velocity 
components  are  given  in  terms  of  <J>  as 

u  «  |i    v  =  |i  (4b) 

dX  '       dy 

2).   Boundary  Conditions 

The  boundary  condition  on  the  bottom  is  that  of  a 
non-porous  wall,  i.e., 

§±  =  0   at  y  =  -h  (5) 

16 


Two  boundary  conditions  are  encountered  at  the  water 
surface,  one  dynamic  and  one  kinematic. 

The  dynamic  boundary  condition  is  obtained  from 
Bernoulli's  Equation  by  setting  the  pressure  equal  to  zero. 
That  is: 

u2  +  v2  +  2||-  -2g(n  -  y)   at  y  =  n        (6) 

where  IT  is  the  total  energy  head,  a  constant. 

Assuming  the  wave  travels  without  change  in  form,  it  is 
then  possible  to  choose  a  reference  system  moving  in  the 
positive  x  direction  with  wave  celerity  C.   This  makes  the 
fluid  motion  steady  with  respect  to  the  moving  reference 
system.   The  Bernoulli  Equation  in  such  a  reference  system 
is: 

(u  -  C)2  +  v2  =  2g(n  -  y)    at  y  =  n        (7) 

The  kinetic  boundary  condition  imposes  the  condition  that 
no  fluid  be  transported  across  the  free  surface.   Or,  in 
other  words,  the  free  surface  must  be  a  streamline.   This  may 
be  expressed  mathematically  in  the  form: 

Equations  (4a),  (5),  (7),  and  (8)  completely  specify  the 

boundary  value  problem. 

There  are,  historically,  several  approximate  solutions 

to  this  problem.   The  solution  to  be  considered  further  in 

this  analysis  is  that  commonly  known  as  Stokes  fifth-order 

wave  theory. 
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D.   STOKES  FIFTH  ORDER  SOLUTION 

Stokes  proposed  [9]  a  power  series  type  of  solution  to 
the  above  problem  and  proved  the  validity  of  the  second  order 
solution  based  on  this  power  series.   Skjelbreia  and 
Hendrickson  [9]  have  extended  this  theory  to  fifth-order, 
yielding  a  solution  to  the  boundary  value  problem  of  the 
form: 

<J>  =  ^-{A  cosh(BS)sin9  +  A2cosh(2BS)sin(2G ) 

+  A-,cosh(33S)sin(39)  +  Ailcosh(4BS)sin(^6) 
+  Accosh(56S)sin(56)}  (9a) 

where  A,  through  A,-  are  as  given  in  Appendix  A. 

This  equation  may  be  written  in  dimensionless  form  as: 

<j>  =  C{A..cosh($S)sine  +  A0cosh(2$S)sin(26 )  +  .  .  .  } 

1  2  (9b) 

where:  <j>   =   0(f>/~\/gh 

C   =   C/Vih 
Differentiation  of  equation    (9a)   with  respect   to  time 


yields 


||-  =  C2{A  cosh(6S)cosG  +  2A2cosh(23S)cos (20 ) 

+   3A3cosh(36S)cos(36)    +   4Aijcosh(4BS)cos(4e) 
+   5A5cosh(5BS)cos(56)  }  (10a) 
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or ,  in  dimensionless  form 

in*  c2 

+  ...}  (10b) 


C  {A1cosh(3S)cos(9)  +  2A2cosh(23S)cos(26) 


where     at  "  at  ^^h" 

The  series  form  of  the  wave  profile  is  assumed  to  be: 

y  =  jkxcos(e)  +  B2cos(2e)  +  B_cos(39) 

+  B^cosC^e)  +  B5cos(56)}  (11a) 

where  B?  through  B,-  are  as  given  in  Appendix  A  and  X  is  an 
arbitrary  constant. 

Equation  (11a)  can  be  expressed  in  dimensionless  form 
as 

y  =  -{X  cos(e)  +  B~cos(2e)  +  ...}  (lib) 

d  * 

where     y  =  y/d 

d  =   27rh/L 
The   series    form  for  II   is   assumed  to  be: 

n  =  |u2c3  +  AUCj,}  (12) 

where  C_  and  Cj,  are  as  given  in  Appendix  A. 

The  wave  celerity  expressed  in  series  form  is 

C  =  £  C2{1  +  X2C,  +  X4C0}  (13a) 

p    O  J.         c. 
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or,  in  dimensionless  form 

C  =  -  C2{1  +  X2Cn  +  X4Cp}  (13b) 

d  ° 

where  C.  and  C?  are  as  given  in  Appendix  A. 

Using  equations  (4b)  and  (9b)  we  find  the  horizontal 

particle  velocity  to  be: 

u  =  (|4)  =  C{A1cosh(BS)cos(6)  +  2A0cosh(2BS)cos(2G) 

d  X  J.  c. 

+   3A3cosh(3$S)cos(3e)   +   4A1|cosh(4BS)cos(4e) 
+   5A5cosh(56S)cos(56)}  (14) 

where  u  =  u/\/gh 

Also,    from  equations    (4b)    and    (9a): 

v  =    (|^)    =   C{Ansinh(BS)cos(0)    +   2A0sinh(2BS)cos(26 ) 
dy  i  c. 

+   3A3sinh    (36S)cos(36)    +   4A1|sinh(4BS)cos  (4e) 
+   5A5sinh(5BS)cos(56)}  (15) 

where  v  =   v/vgh 

Differentiation  of  equations  (14)  and  (15)  with  respect 
to  time  yields  the  particle  accelerations 

ax  =  (|^)  =  dC2{A1sinh(BS)sin(6)  +  2A2sinh(2BS)sin(29) 

+   3A3sinh(3BS)sin(36)    +   4A|jSinh(4BS)sin(4e) 

+   5Acsinh(5BS)sin(56)  (16) 

z> 

where  ax  =  ax/g 
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and 


a   =  dC2{A1sinh(3S)cos(9)  +  2A2sinh(2eS)cos(26) 
+  ...}  (17) 


where     a  =  a  /g 

y   y 

The  potential  function  and  it's  derivatives,  as  well  as 
the  wave  profile  are  now  known  in  terms  of  the  unknown 
quantities  X  and  d.   Specifying  equations  for  these  parameters 
are  found  by  developing  equations  for  the  peak-to-trough 
wave  height  and  the  wave  period.   That  is,  the  dimensionless 
wave  height  is  defined  as: 

Hd  =  ?e=0  -  J6=u  (18) 

where     Hd  =  H/h 

Substitution  of  equation  (13b)  into  equation  (18 ) 
yields : 

Hd  =  2{X  +  X3B33  +  X5(B35  +  B55)}  (19) 

where  B~~,  B,-,  and  B,-,-  are  as  given  in  Appendix  A. 

Now,  the  wave  period  and  the  wave  celerity  are  related 
by  the  relationship: 

C  =  L/T  (20) 


So,  equation  (15a)  becomes 

C2 
(L/T)2  =  g-^{l  +  X2C,  +  X4CP}  (21) 

6         L  d 
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or  , 

h/gT2    =        d0      C2(l   +    X2C,    +    xV0}  (22) 

2|IX  °  x 

2 

The  parameter  gT  /h  is  referred  to  herein  as  the  period 

parameter. 

Equations  (19)  and  (22),  if  solved  simultaneously,  yield 

2 
values  for  X  and  d  in  terms  of  Hd  and  gT  /h,  quantities  which 

are  incident  wave  parameters.   Thus  the  potential  function, 

particle  velocities  and  accelerations  are  shown  to  be  functions 

of  the  dimensionless  parameters  obtained  by  dimensional 

analysis  and  given  by  equation  (3).   (Note  that  Hd  =  2H/h) 

The  next  step  in  the  analysis  is  to  express  the  forces 

acting  on  the  cylinder  as  functions  of  the  potential  function, 

particle  velocities  and  accelerations. 

E.   FORCES  ON  THE  CYLINDER 

The  forces  acting  on  the  cylinder  are  described  in  terms 
of  their  horizontal  and  vertical  components. 

1).   Horizontal  Component 

The  horizontal  component  of  force  acting  on  the 
cylinder  is  expressed  in  the  form  of  the  so-called  Morrison 
Equation  as 

Q 

Fx  =  -r(p2a£u2)  +  (1.0  +  CM)7ra2*ax  (23) 

where     C~  =  drag  coefficient 

C..  =  added  mass  coefficient 
M 
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For  cases  where  fluid  particle  motion  is  small  in 
comparison  with  the  cylinder  diameter,  flow  separation  does 
not  occur  and  the  contribution  of  drag  to  the  total  force 
may  be  disregarded.   In  this  case,  equation  (23)  may  be 
rewritten 

Px  =  (1.0  +  CM)ax  (24) 

2 
where     F  =  F  /pga  I. 

The  added  mass  coefficient  for  a  circular  cylinder  in 

contact  with  a  rigid  wall  was  given  in  closed  form  by 

Garrison  [13]  as  CM  =  2.29. 

2).   Vertical  Force  Component 

As  mentioned  previously,  application  of  the  Morrison 

Equation  to  the  vertical  component  of  force  yields  a  result 

which  is  clearly  erroneous.   The  vertical  force  may  be 

considered  to  arise  from  two  sources:  an  inertial  force 

and  a  lift  force. 

a).   Lift  Force 

The  lift  force  may  be  expressed  as: 

FL  =  l/2(2aAp)CLu|u|  (25) 

where     CT  =  lift  coefficient. 

Li 

or,  in  dimensionless  form: 

FL  =  hCLu|u|  (26) 

2 

where     F,  =  F,./pga  I. 
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Dalton  and  Helf stein  indicate  [7]  that  C  =4.49 
for  unseparated  flow  past  a  circular  cylinder  in  contact 
with  a  rigid  wall. 

b).   Inertial  Force 

The  pressure  distribution  around  the  surface  of 
the  cylinder  existing  if  the  cylinder  is  not  present  results 
in  a  net  vertical  force  (known  as  a  Froude-Krilov  force). 

Consider  the  resultant  differential  force  dP. 
acting  on  a  differential  area  of  the  cylinder  dA  as  shown 
in  Figure  2.   The  above  equation,  written  in  dimensionless 
form,  is: 

F  =   /    P(a,Y)sinydY  (30) 

L   0 

The  pressure  is  obtained  from  the  Bernoulli  Equation: 

P  =  -p{|(u2  +  v2)  +  ||}  (3D 

or,  in  dimensionless  form: 

P  =  -h[|(u2  +  v2)  +  (||)]  (32) 

where     P  =  P/pga 
Therefore 

i±   =  -h  /[|(u2  +  v2)  +  (f£)]sinYdY         (33) 


And  the  total  vertical  force  may  be  written: 

2n 

'      I 
0        c  ou  (34) 


P_   =   h{CTu|u|    -   /      [i(u2   +   v2)    +    (|i)]sinYdY} 
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INERTIAL    FORCE    GEOMETRY 
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F.   DISCUSSION  OF  COMPUTER  PROGRAM 

The  equations  developed  above  were  programmed  for  a 

digital  computer  as  given  in  Appendix  D.   The  basic  inputs 

~  2 

to  the  program  are  H,  h  and  gT  /h.   The  program  yields 

particle  velocity,  acceleration  and  computes  the  horizontal 

and  vertical  forces  on  the  cylinder  at  various  wave  angles. 

The  STAKES5  subroutine  calculates  the  coefficients  for 
the  Stokes  fifth  order  equations.   Simultaneous  solution 
to  equations  (19)  and  (22)  are  obtained  by  assuming  an  initial 
value  of  X  and  d  based  on  linear  wave  theory  and  iterating 
to  the  actual  value  using  the  Newton-Raphson  technique. 

The  main  program  generates  tables  of  values  and  in-line 
graphs  of  the  wave  profile,  particle  velocity  vs.  depth, 
vertical  and  horizontal  force  vs.  phase  angle. 
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III.   DESCRIPTION  OF  APPARATUS  AND 
EXPERIMENTAL  PROCEDURE 

A.   WAVE  CHANNEL 

Wave  forces  on  the  cylinder  were  measured  experimentally 
in  a  wave  channel.   The  basic  wave  channel,  as  shown  in 
Figure  3,  was  fifteen  inches  wide  and  consisted  of  two  sec- 
tions: a  transition  section  and  a  shallow  water  section. 
The  overall  channel  length  was  87  feet. 

Three-quarter  inch  exterior  plywood  sheets  were  used  for 
the  channel  sides.   Vertical  rigidity  was  provided  by  bracing 
the  plywood  sheets  with  2x4  fs  at  various  intervals.   The 
vertical  members  were  placed  four  feet  on  centers  in  the 
shallow  water  section  and  from  four  feet  to  eighteen  inches 
in  the  transition  section.   Additional  rigidity  was  provided 
in  the  transition  section  by  truss  arrangements,  as  shown 
in  Figure  4.   Two  2x4 's  were  attached  to  the  top  of  the 
vertical  studs  and  extended  the  length  of  the  channel. 

A  double  floor  constructed  of  two  plywood  sheets  with 
two  half-inch  separators  was  used.   The  sides  of  the  channel 
were  bolted  together,  through  the  spacing  of  the  double 
bottom,  with  three-eights  inch  threaded  rod. 

All  wooden  sections  exposed  to  water  were  water-proofed 
and  painted  with  three  coats  of  Morewear  Vitri-Glaz  1320A. 
Dow  Corning  Sealant  780  was  applied  to  all  seams  and  joints. 

A  paddle  type  wave  generator  was  used  to  generate  waves 
in  the  channel.   An  aluminum  plate  was  hinged  at  the  tank 
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Fig. 3      WAVE  CHANNEL 
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Fig.4 
TRANSITION   SECTION    TRUSS  ARRG'T. 
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floor  and  attached  by  a  pin  connection  at  the  top  to  a  driving 
rod.   Wave  motion  behind  the  paddle  was  damped  by  use  of 
baffle  plates.   A  two  horse  power  variable  speed  drive,  with 
an  output  speed  range  of  twenty  to  one  hundred  eighty  revolu- 
tions per  minute,  was  mounted  on  top  of  the  wave  channel  and 
fitted  with  a  six  and  three  quarter  inch  radius  face  plate. 
Ways  were  attached  to  the  plate  to  adjust  the  driving  rod 
eccentricity  between  zero  and  eight  inches. 

The  transition  section  contained  a  ramp  with  a  slope  of 
1:20  which  shoaled  the  generated  wave  with  a  change  in  depth 
from  four  to  two  feet.   The  overall  length  of  the  transition 
section  was  42  feet. 

The  shallow  water  section  had  a  depth  of  two  feet  and 
an  overall  length  of  42  feet.   The  test  module  was  located 
twenty  four  feet  from  the  transition,  allowing  the  waves  to 
reach  a  fully  developed  state  following  the  transition  before 
reaching  the  test  section. 

A  variable  slope  dissipater  beach  was  located  at  the  end 
of  the  test  section  in  order  to  dissipate  the  wave  motion. 
The  beach  was  constructed  of  metal  shavings  held  between  two 
pieces  of  perforated  stainless  steel  sheet  metal,  which  were 
separated  by  two  inch  wooden  spacers.   A  solid  sheet  of  one- 
eighth  inch  aluminum  was  attached  to  the  bottom  of  the  beach 
by  three  and  five  eighths  inch  separators.   One  by  one  inch 
aluminum  angles  were  attached  to  the  exposed  surface  of  the 
beach,  parallel  to  the  wave  fronts.   Maximum  slope  of  the 
beach  was  1:7. 
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B.   TEST  MODULE 

The  test  module  consisted  of  the  circular  cylinder  with 
its   own  floor  and  walls,  constructed  as  an  integral  system 
to  provide  easy  removal  of  the  model  from  the  tank.   The  wave 
channel  sides  and  floor  were  recessed  to  provide  a  smooth 
transition  from  channel  to  test  module.   The  channel  walls 
were  made  of  plexiglass  in  the  test  module  region  to  allow 
maximum  visibility  of  the  model.   The  test  module  is  shown 
in  Figure  5. 

A  four  inch  diameter  plexiglass  cylinder  was  supported 
between  the  walls  of  the  test  module  by  use  of  cantilever 
beams  and  adjusted  to  approximately  one-sixteenth  inch  above 
the  floor.   A  thin  flexible  plastic  barrier  was  installed  in 
this  gap  to  prevent  any  water  motion  under  the  cylinder  due 
to  wave  motion.   The  barrier  was  held  in  place  by  '0'  ring 
material  pressed  into  0.007  inch  slots  in  the  cylinder  and 
test  module  floor. 

The  cantilever  beams  were  used  to  support  the  cylinder, 
as  shown  in  Figures  6  and  7.   Bulkheads  were  fixed  in  the 
cylinder  at  approximately  two  inches  from  each  end,  and  the 
fixed  ends  of  the  cantilever  beams  bolted  on  these  bulkheads. 
The  free  ends  of  the  two  beams  protruded  slightly  beyond  the 
end  of  the  cylinder  and  into  the  plexiglass  test  module  walls 
where  they  were  supported  by  small  self-aligning  ball  bearings 
pressed  into  the  plexiglass.   Both  beams  were  fitted  with 
strain  gauges  and  waterproofed  using  BLH  Barrier  '  C,  as 
shown  in  Figure  8. 
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Fig.  5   TEST  MODULE 
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Fig.  6   HORIZONTAL  FORCE  CANTILEVER 
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Pig.  7   VERTICAL  FORCE  CANTILEVER 
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Fig.    8      FORCE    CANTILEVER 
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One  beam,  positioned  with  its  largest  cross  sectional 
dimension  vertical,  was  used  to  measure  horizontal  force  and 
the  other,  with  its  largest  cross  section  positioned  hori- 
zontally, was  used  to  measure  vertical  force.   Calibration 
tests  indicated  that  cross  coupling  was  negligible. 

Dimensions  of  the  beams  were  chosen  to  provide  sufficient 
flexibility  to  allow  measurement  of  forces  using  strain  gauges 
and  sufficient  stiffness  to  provide  a  natural  frequency  that 
was  large  in  comparison  with  the  excitation  frequency.   In 
addition,  the  minimum  width  of  the  beams  was  constrained  by 
the  width  necessary  to  mount  the  strain  gauges.   A  reasonable 
trade-off  was  arrived  at  by  making  use  of  the  maximum  sensi- 
tivity of  the  amplifier/recorder  and  determining  a  maximum 
allowable  value  for  the  cantilever  beam  length. 

To  prevent  rotation  of  the  cylinder,  an  offset  arm  was 
used  to  connect  the  horizontal  force  cantilever  to  the  test 
module.   A  Farber  Bearing,  AMS  1K7,  was  used  as  a  wheel  at 
the  test  module  end  of  the  arm  to  reduce  friction. 

C.   WAVE  HEIGHT  PROBE 

A  capacitance  type  probe  was  chosen  to  measure  wave 
heights  in  this  study.   The  probe  essentially  consists  of 
a  single  insulated  wire  which  acted  as  one  plate  of  a  capac- 
itor, the  water  acting  as  the  second  plate.   The  capacitance 
varied  as  the  water  level  rose  and  fell.   A  three-eighths 
inch  diameter  acrylic  rod  was  attached  to  the  bottom  of  the 
foil-like  support  member  and  the  sensing  wire  was  attached 
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at  the  bottom  of  the  acrylic  rod  and  insulated  there.   The 
tip  of  the  acrylic  rod  was  removable  to  provide  easy  changing 
of  the  sensing  wire.   The  wire  was  connected  at  the  top  to 
a  cable  connector  mounted  on  a  nylon  block.   Number  30  A.W.G. 
wire  with  poly  thermaleze  insulation  was  used. 

The  schematic  of  the  electronic  circuit  utilized  with 
this  probe  is  shown  in  Figure  9.   In  this  circuit,  the 
square  wave  output  of  a  multi-vibrator  is  used  to  drive  a 
capacitance  bridge,  one  leg  of  which  is  varied  by  the  wave 
height  probe.   The  resulting  AC  error  signal  is  converted  to 
a  D.C.  output  signal  by  the  full  wave  rectifier.   The  output 
signal  was  found  to  be  insensitive  to  small  multivibrator 
frequency  changes  and  linear  through  the  range  of  the  wave 
height  probe.   A  model  LD  5211-13  Brush  D.C.  amplifier  was 
used  with  the  probe  circuit. 

D.   TEST  PROCEDURE 

Calibration  of  the  probe  was  conducted  with  the  tank 
filled  to  the  desired  level  and  the  test  module  and  wave 
height  probe  in  position.   With  the  probe  immersed  to  a  depth 
of  nine  inches,  the  amplifier  gain  was  set  at  desired  sensi- 
tivity and  the  recorder  pen  bias  adjusted  to  position  the 
recorder  to  mid-scale.   The  probe  immersion  was  varied  over 
a  seven  inch  range  in  half  inch  increments,  using  a  traverse 
mechanism.   This  initial  calibration  was  checked  at  several 
points  before  and  after  each  set  of  runs. 
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Fig.  9    WAVE    HEIGHT    PROBE  CIRCUIT 
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The  force  cantilevers  were  calibrated  by  loading  the 
cylinder  with  a  series  of  weights,  using  a  pulley  arrange- 
ment shown  in  Figure  10.   The  load  was  transmitted  to  the 
cylinder  through  a  series  of  tapped  holes  located  around  the 
circumference  at  positions  horizontally  fore  and  aft,  and 
vertical.   The  Brush  recorder  output  was  calibrated  by  loading 
the  cylinder  in  half-pound  increments.   This  calibration 
was  also  checked  prior  to  and  after  each  series  of  runs. 

After  calibration  was  completed,  the  wave  generator  speed 
was  set  for  the  desired  wave  length  and  a  data  run  commenced. 
During  each  run,  wave  generator  eccentricity  was  varied  from 
one  inch  to  eight  inches,  or  to  a  point  where  the  waves 
broke  in  the  channel. 
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FIGURE  1Q. CANTILEVER  CALIBRATION  ARRANGEMENT 
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IV.   PRESENTATION  OF  RESULTS  AND  CONCLUSIONS 

As  discussed  in  the  theoretical  analysis,  dimensionless 
wave  force  coefficients  can  be  represented  as  functions  of 
the  dimensionless  water  depth  (h/a),  the  dimensionless  wave 
height  (H/2a),  and  dimensionless  period  (gT  /h)  parameters. 
A  test  program  was  developed  which,  for  a  given  run,  held 
water  depth  and  period  constant  while  varying  wave  height, 
in  order  to  determine  the  effect  of  the  dimensionless  param- 
eters on  the  horizontal  and  vertical  force  coefficients. 
Four  series  of  experimental  runs  were  conducted  at  relative 
water  depths  of  9.0,  7.0,  5.5,  4.0. 

The  dimensionless  wave  period  parameter  was  varied  for 
each  series  over  a  range  of  100  to  the  maximum  value  attain- 
able with  the  experimental  apparatus.   The  upper  limit  in 
wave  height  was  usually  provided  by  breaking. 

A  series  of  representative  traces  are  shown  in  Figures 
11,  12,  13  and  14  due  to  the  difficulty  in  presenting  all 
the  traces.   These  traces  represent  experimental  data  for 
three  dimensionless  periods  from  each  of  the  four  dimension- 
less water  depths.   The  uppermost  trace  in  each  run  repre- 
sents wave  height,  the  middle  trace  represents  vertical 
force  and  the  bottom  trace  represents  horizontal  force.   The 
experimental  values  of  the  force  coefficients  for  each 
experimental  case  are  given  In  Appendix  B. 

Using  the  dimensionless  period,  water  depth  and  waveheight 
as  inputs,  the  computer  program  was  used  to  generate 

Ml 


o 

II 

a 


w 

(D 
O 

a 
u 

EH 

<l) 
o 

O 

Ph 

I 

-p 

bO 
•H 
<D 

0) 

> 


C5 
M 


42 


t 

^, 

[It 

I 

\ 

i\i  ' 

i  \ 

i  j 

1 

X 

1 

i  r 

p 

___ 

1  i 

U^- 

i 

'     Is- 

H — 

i^- 

* 

2%Z 

*^ 

*~     ~ i 

H 

V 

V 

s 

\ 

/ 

T 

_^ 

' 

/ 

' 

:i 

\ 

J 

I 

l 

jr 

T 

_i 

1  f 

I     3 

-  4X 

W 

i 

i 
s 

ZST 


*~"~v~ — — 


K'H-]-t± 


II 

cd 
xi 


w 

0) 

o 
a 
u 

Eh 

d) 
o 

u 

O 

[*« 

I 
-P 

w 

•H 
<D 
K 

<D 
> 

a 


CO 


C5 
M 


43 


} 

*  /  _. 

, 

in       \ 
o         > 

in     v 

H 

J?i   Y,  1 

■*       /     1 

< 

■  <S      j 

44 


. 

n     -y 

!  .  y. 

--:v~- 

»_ 

____ 

\ 

I 

^    /] 

u>  • 

& 

K 

3 

a: 


45 


corresponding  theoretical  values  of  the  force  coefficients 
for  each  experimental  case.   This  provided  a  comparison 
between  experimental  and  theoretical  values  for  each  run. 
Plots  of  experimental  and  theoretical  values  for  maximum 
horizontal  and  vertical  force  coefficients  as  a  function  of 
wave  height  for  each  period  and  water  depth  studied  are 
included  in  Appendix  C. 

Experimental  results  and  theoretical  computations 
indicate  that,  at  the  lower  values  of  dimensionless  period 
and  wave  height  studied,  the  wave  may  be  characterized  as  a 
linear  wave.   The  horizontal  force  was  found  to  vary  sinu- 
soidally,  and  increased  linearly  with  waveheight.   The 
vertical  force  varied  approximately  as  a  negative  cosine 
function  in  this  region,  having  its  maximum  value  in  the  wave 
trough.   The  magnitude  of  the  vertical  force  coefficient 
was  found  to  vary  approximately  as  the  waveheight  squared. 
These  characteristics  are  shown  most  clearly  in  runs  Number 
44  through  47,  71  through  74  and  93  through  95  in  Appendix  C. 

As  the  relative  period  and  wave  height  increased,  the 
wave  displayed  the  sharper  peaks  and  longer,  shallower  troughs 
that  are  characteristic  of  non-linear  gravity  waves.   The 
horizontal  force,  proportional  to  the  fluid  particle  acceler- 
ation, tended  to  reach  maximum  values  closer  to  the  wave 
peaks,  but  still  remained  symmetrical,  with  respect  to  the 
wave  crest,  as  indicated  by  the  theoretical  traces.   The 
theoretical  analysis  predicted  a  phase  lag  between  the  wave 

o 

crest  and  the  maximum  force  as  small  as  24   at  the  higher 

wave  heights. 
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The  vertical  force  component  began  to  display  positive 
values  occurring  at  twice  the  frequency  of  the  wave  peaks, 
as  shown  in  the  experimental  traces.  Figures  11,  12,  13,  1*1 . 
This  is  explained  by  the  increasing  importance  of  the  lift 
component  of  vertical  force,  which,  theoretically,  is  similar 
to  a  cosine  squared  function.   At  extremely  long  periods,  the 
lift  force  is  clearly  dominant,  and  the  maximum  force  no 
longer  occurs  in  the  middle  of  the  wave  trough.   The  maximum 
vertical  force  occurs  when  the  cylinder  is  under  the  wave 
crests  where  the  velocity  is  greatest. 

At  large  periods  and  wave  heights,  a  point  is  reached 
where  flow  separation  effects  are  no  longer  negligible.   This 
is  seen  most  clearly  in  the  vertical  force  component,  whose 
value  becomes  suddenly  lower  than  the  theory  predicts,  though 
still  varying  as  the  waveheight  squared.   This  sudden  drop 
in  lift  is  apparently  caused  by  a  corresponding  decrease  in 
the  lift  coefficient  resulting  from  flow  separation.   This 
phenomenon  is  most  clearly  evident  in  the  graphical  presentation 
of  runs  51  through  56,  and  71  through  73  in  Appendix  C. 

1).   Horizontal  Forces 

As  previously  mentioned,  the  horizontal  forces 
tended  to  be  nearly  sinusoidal  in  nature  with  a  maximum  value 
which  is  quite  linear  with  wave  height.   It  is  therefore 

possible  to  characterize  the  horizontal  force  coefficient 

2 

by  it's  maximum  value  (F    /pga  £).   The  variation  in  the 

xmax 
horizontal  force  coefficient  with  dimensionless  wave  height 

is  shown  graphically  in  Appendix  C. 
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The  experimental  values  of  horizontal  force  are  seen 
to  vary  linearly  with  wave  height  for  relative  period  param- 
eter values  from  100  to  120.   For  relative  period  values 
greater  than  120,  the  force  coefficient  varies  in  an  increas- 
ingly non-linear  fashion.   This  is  consistent  with  the  theory, 
as  shown  in  Figure  15.   However,  for  the  range  of  period 
parameters  studied,  the  deviation  from  linearity  is  relatively 
small.   It  is  therefore  possible  to  display  the  experimental 
data  in  terms  of  the  slope  of  the  force  coefficient  versus 
wave  height  plots,  as  shown  in  Figures  16  and  17. 

The  theory  predicts  the  maximum  values  of  horizontal 
force  coefficients  prior  to  the  wave  heights  where  separation 
effects  become  apparent.  The  discrepancy  between  theoretical 
and  experimental  values  is  due  to  variation  of  the  added  mass 
coefficient  in  a  separated  flow  as  well  as  a  net  drag  force 
which  is  no  longer  zero. 
2).   Vertical  Forces 

As  indicated  previously,  the  vertical  force  displays 
two  regimes;  one  inertia  dominated  and  one  lift  dominated. 

The  maximum  value  of  the  vertical  force  coefficient 

2 
(F    /pga  I)    is  seen  to  vary  as  the  wave  height  squared 

^max 

throughout  the  range  of  data  taken.   This  is  shown  graphically 

in  Appendix  C.   For  higher  values  of  relative  wave  height, 

the  magnitude  of  the  vertical  force  coefficients  is  less  than 

predicted  by  the  theory,  due  to  separation  effects,  as 

explained  above. 
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Fig.  15     HORIZ.  WAVE    FORCE 
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B.   CONCLUSIONS 

From  the  experimental  and  theoretical  results,  the 
following  conclusions  are  considered  warranted: 

1.  The  horizontal  forces  on  the  cylinder,  resulting 
from  the  incident  gravity  waves,  increased  linearly  with 
wave  height  for  relative  period  values  less  than  120.   For 
relative  period  values  greater  than  120,  the  force  coeffi- 
cient was  found  to  vary  with  wave  height  in  an  increasingly 
non-linear  fashion. 

2.  The  vertical  forces  displayed  two  regimes;  one  inertia 
dominated  at  lower  periods  and  wave  heights,  and  one  lift 
dominated  at  higher  values  of  these  parameters. 

3.  For  high  period  parameter  values,  separation  effects 
can  be  expected  for  relative  wave  height  values  above  0.6, 
resulting  in  a  sudden  drop  in  experimental  force  coefficient 
values  from  those  predicted  by  the  theory. 

h .   The  theory  used  in  this  study  provides  excellent 

2 
agreement  with  experimental  values  for  relative  period  (gT  /h) 

values  from  100  to  200,  and  prior  to  the  onset  of  flow 

separation.  . 
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APPENDIX  A 
COEFFICIENTS  FOR  STOKES  FIFTH  ORDER  WAVE  THEORY 
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APPENDIX  B 
EXPERIMENTAL  DATA 
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APPENDIX  C.   WAVE  HEIGHT-FORCE  PLOTS 
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